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We calculate a current and its fluctuation in a two-state stochastic system under a periodic 
perturbation. The system could be interpreted as a channel on a cell surface or a single Michaelis- 
Menten catalyzing enzyme. It has been shown that the periodic perturbation induces so-called pump 
current, and the pump current and its fluctuation are calculated with the aid of the geometrical 
phase interpretation. We give a simple calculation recipe for the statistics of the current, especially 
in a non-adiabatic case. The calculation scheme is based on the non-adiabatic geometrical phase 
interpretation. Using the Floquet theory, the total current and its fluctuation are calculated, and it 
is revealed that the average of the current shows a stochastic-resonance-like behavior. In contrast, 
the fluctuation of the current does not show such behavior. 
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I. INTRODUCTION 

Recently, it has become possible to perform single- 
molecule experiments. [l| In such single-molecule ex- 
periments or small chemical systems such as cells, it is 
expected that 'fluctuation' plays an important role be- 
cause mean values and fluctuations (deviations) of ob- 
servables are in the same order. 0, H| It also becomes 
possible to measure flux distributions experimentally in 
such small systems. [|[ When the fluctuation is large, the 
traditional rate equation approach is not adequate for the 
study of stochastic systems because it treats only average 
values in large size limit. Hence, the stochastic behavior 
of chemical reaction systems should be treated by a mas- 
ter equation approach. [f| Q Although it is difficult to 
obtain exact solutions of master equations in general, it 
is important to evaluate (even approximately) not only 
the mean values, but also the fluctuations. 

There is a phenomenon, known as a pumping, in which 
a system under a periodic perturbation causes a finite 
flux in a preferred direction. For example, a classical 
Michaelis-Menten (MM) enzymatic mechanism under a 
periodic perturbation has been studied experimentally 
and theoretically. It has been revealed that the periodic 
perturbation activates a pumping mode in the MM sys- 
tem. 0, H, Q That is, the periodic perturbation causes 
a current which is not explained by a simple average of 
those in the strict static cases. The concept of the pump 
current is also related to molecular motors and Brownian 
ratchets. [TSUI- 

The average current can be calculated by various meth- 
ods. [1, [H, [H, [3| In addition, it has been indicated that 
a calculation of the net current in the Brownian ratch- 
ets shows a similar mathematical structure to that of the 
geometric phase in quantum mechanics. [l5[ Recently, 
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an explicit connection between the statistics of the cur- 
rent flow and the Berry phase interpretation has been 
proposed. [TtJ The total current is divided into the 
so-called 'classical' current and 'pump' current. Sinit- 
syn and Nemenman have shown that the dynamical and 
the Berry phase correspond to the classical and pump 
current, respectively. [16j | One of the other remarkable 
progresses is that it becomes possible to calculate not 
only the average current, but also the fluctuation about 
the current for an 'adiabatic' cases. Here, the 'adiabatic' 
means that the oscillation of the periodic perturbation is 
very slow. 

Although the average current has already been ob- 
served experimentally, [7[ any experiment for the fluctu- 
ation has not been performed. However, recent progress 
of experimental techniques Q would make it possible to 
observe the flux distribution for the pumping phenom- 
ena directly. Detailed information for non-equilibrium 
systems is useful to construct the non-equilibrium sta- 
tistical mechanics and non-equilibrium thermodynamics. 
18] In order to investigate such non-equilibrium systems 
experimentally, theoretical predictions would become a 
guideline for experimentalists. 

In the present paper, we give a calculation scheme for 
the statistics of the current, especially in 'non-adiabatic' 
cases. The fluctuation about the average current in the 
non-adiabatic cases has not been calculated yet. In ad- 
dition, the non-adiabatic region is of great interest be- 
cause a stochastic resonance-like phenomenon has been 
observed experimentally in such non-adiabatic region. Q 
One of our main results is that there is no stochastic 
resonance-like behavior for the fluctuation about the av- 
erage current. In order to calculate the statistics, we 
develop a calculation scheme based on the non-adiabatic 
geometrical phase (Aharonov-Anandan phase [lfl [20|). 
We have already proposed a method for direct evaluation 
of the Aharonov-Anandan phase by means of a pertur- 
bation calculation for eigenvectors of a Floquet Hamilto- 
nian. 21] Here, we show that a perturbation calculation 
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for eigenvalues of a Floquet Hamiltonian gives the total 
statistics directly; the calculation scheme is simpler than 
that in Ref. 21. 

The present paper is organized as follows. In Sec. II, 
we explain the stochastic model for the pumping phe- 
nomenon, which has been introduced by Sinitsyn and 
Nemenman. fl6| In Sec. Ill, we review a geometrical 
phase interpretation for the pumping phenomenon. We 
develop an analytical treatment for the calculation of the 
total current with the aid of the Floquet theory and the 
perturbation calculation in Sec. IV. Section V gives con- 
cluding remarks. 



II. STOCHASTIC PUMPS 

We mainly focus on a simple stochastic system with 
two absorbing states, which has been introduced in Ref. 
16. The classical stochastic system is described as fol- 
lows: 



[L] t=> [C] «=> [R]. 



(1) 



The system consists of three parts. Two absorbing states 
are denoted as [L] and [R]. These absorbing states may 
be interpreted in many ways: they correspond to sub- 
strate and product in a MM enzymatic reaction, or cel- 
lular compartments, and so on. These absorbing states 
exchange molecules or particles via an intermediate con- 
tainer [C]. The container [C] can contain either zero or 
one particle in it. When the container is filled with one 
particle, the particle can escape from the container by 
jumping into one of the two absorbing states, [L] or [R]. 
On the contrary, when the container is empty, either of 
the absorbing states can emit a new particle into the 
container. For simplicity, we here consider the following 
kinetic rates: 



(i) [L] -> [C] : ki = ci + i?cos(wi), 

(ii) [L]^[C]: fc_i=c_i, 

(iii) [C] [R] : k 2 = c 2 , 

(iv) [C] «- [R] : fc_ 2 = c_2 + Rsm{ojt), 



(2) 



where c±\ and c±2 are positive real numbers, and R is 
the amplitude of the perturbative oscillation. The above 
kinetic rates indicate that only two kinetic rates (ki and 
fe_2) oscillate with time at a frequency lu. 

Note that the current j at static cases (u> = 0), in 
which all kinetic rates are constant, is easily calculated 
from [l(| 



J = 



K + 



K 



K 



K± = C±lC±2. 



(3) 



{m} 



We here describe the current j as classical current. One 
might think that the time average of the current over a 
cyclic perturbation is simply given by the time average 
of the classical current j, but it is not true. For example, 
when c\ = c_i = c 2 = c_ 2 and R ^ 0, the time average 



of j gives zero. However, even in the simple case, a net 
current is actually observed. Such additional current is 
called a pump current. The problem considered here is to 
calculate the total current which includes both classical 
and pump currents. 



III. PHASE INTERPRETATION FOR THE 
CURRENT STATISTICS 

A. Generating function for the current 

Our main goal is to calculate the net current between 
[C] and [R] in the steady state. In order to calculate 
the current, we use a method similar to the full count- 
ing statistics. [22] Let P n be the probability to have n 
net transitions from [C] into [R] during time T, where 
T = 2ir/oj is the period of the rate oscillations. The 
probabilities of a filled and empty state of [C] are de- 
noted by Pf and P c respectively, and the state of the 
system is defined by 



p(*) 



Pa 
Pi 



(4) 



Due to the normalization condition, P c + P{ = 1. By 
discussions in Refs. 16 and 23, the characteristic function 
of P n is given by 



Z( X ) = e s ^ = P n 

S— — QO 



(5) 



where 1 is the unit vector, T stands for the time-ordering 
operator, and 



H( X ,t) 



ki + fc-2 
—ki — k-2& 



-k-i — kze %x 
k-i + k 2 



(6) 



In eq. (O, X is called the counting field, and the 
derivatives of S(x) give cumulants of P„, e.g., (n) = 
— idS(x)/dx\ x =o- Hence, the problem for calculating the 
current is actually the evaluation of the characteristic 
function Z(x) (or S{x))- 



B. Interpretation as a Shrodinger-like equation 

From eq. ©, it is easy to see that the characteristic 
function Z(x) is related to a solution of the following 
differential equation: 



-p(t) = -H( X ,t)p(t), 



(7) 



because the final state p(T) at time T is formally calcu- 
lated from 

p(T) = f (e- fo p( ) = expMx)] P (0). (8) 
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Using the normalization condition of the probability 
l^p(O) = 1, we obtain the identity S(x) — *Mx)- Hence, 

(2) 

the total current Jtotai and its derivative J total are given 
by 



^total 



1 dnix) 
T 



J, 



(2) 



total 



H) 



1 9V(X) 



T d X 2 



(9) 
(10) 



Although the existence of such cyclic state has not been 
mathematically justified in general cases, we checked nu- 
merically that the stochastic system reaches a steady 
cyclic state, starting from arbitrary initial conditions. 
Note that when x = 0, eq. ([7]) simply gives a master 
equation for the time evolution of the system, and hence 



n(x = 0) = 0, 



(11) 



due to the cyclic evolution of the system. 

We here use an analogy between the classical stochas- 
tic system and a quantum mechanical formulation. Re- 
placing the time evolution operator H{x, t) by H = 
—iH(xjt), we obtain the following Shrodinger-like equa- 
tion: [24j 

ij t \4>(t))=Hm))- as) 

The time evolution operator U(t) is constructed as 

m) = u(t)\m), as) 

and the time evolution operator satisfies 

ijU(t) = HU(t). (14) 

When we take the initial state \<fi(0)) as a cyclic state, we 
have 

\ct>(T)) = U(T)\cl>(0))=e iflM \mh (15) 

where fi(x) is a phase factor. Hence, the state vector 
\<f)(t)) is an eigenstate of the time evolution operator U(t), 
and the phase factor /z(x) is related to its eigenvalue. 

From this replacement, the problem of calculation of 
the total current and its fluctuation is replaced by the 
evaluation of the phase factor ^(x) for eq. (fTSf . Note 
that the 'Hamiltonian' H is a non-Hcrmitian operator, 
which is different from the usual quantum mechanics. 



IV. DIRECT EVALUATION OF THE TOTAL 
CURRENT AND ITS FLUCTUATION 

In order to calculate the total statistics, the Floquet 
theory is available. [IEIHj E3 The original Floquet the- 
ory for the geometrical phase was applied to Hermitian 



Hamiltonian case s. 1261 12711 and it was extended to non- 
Hermitian cases. 

mm 

As shown in Ref. 21, the 'pump' current is calculated 
from the eigenvectors of the Floquet states. In contrast, 
we here shows that the 'total' current is calculated from 
the eigenvalues of the Floquet states directly. The re- 
lationship among the total phase, the dynamical phase, 
and the Aharonov-Anandan phase is shown in the Ap- 
pendix. 



A. Usage of the Floquet theory 

From the Floquet theorem, the non-unitary time- 
evolution operator U (t) for a periodic non- Hermitian 
Hamiltonian H with period T is decomposed into the Flo- 
quet product form as U(t) — V(t) exp(iMt). Here, U(t) 
is the unique fundamental matrix satisfying t/(0) = /, 
and the non-unitary matrix V(t) also has a period T. 
In addition, M is a time-independent matrix. By using 
the above Floquet product form, the total phase fi(x) is 
related to the eigenvalues of U (T) via 

U(T)\cf> a (0)) = e 4M7 > Q (0)) = e^W|0 a (O)), (16) 

where we define the initial state \<fi a (0)) as the eigenvec- 
tor of U(T). The index a specifies an eigenvector and 
its corresponding eigenvalue. Equation (|16|) means that 
/j, a (x)/T is the eigenvalues of M. From the above discus- 
sion, it is needed to calculate the time evolution operator 
U(t) or the matrix M and then the evaluation of its eigen- 
values should be done. However, it is difficult to obtain 
them in general. 

In order to evaluate the total phase A* a (x), we here 
introduce another fundamental matrix and its Floquet 
product form F(t) = P(t) exp(iQt), in which Q is as- 
sumed to be a diagonal matrix. [2l|,[23| The fundamen- 
tal matrix F(t) is useful for the Fourier analysis. Since U 
and F are both fundamental matrices of H, there exists 
a constant invertible matrix X with U(t) = F(t)X. (2?| 
Because U(0) = I, we have U(0) = F(0)X = /. Hence, 
we obtain X — F(0) -1 and 

U(t) = F^FiO)- 1 = P(t)P(0)- 1 e 4P(0)QP(orl *. (17) 

Here, we can make identifications V(t) — P(t)P(0) _1 
and M = P(Q)QP(0)^ 1 . Using the above identifications 
and 



M\<p a (0)) 



Mo(x) 
T 



16,(0)), 



we obtain 



Q\a) = - \a) 



(18) 



(19) 



where \a) = P(0) -1 |</> a (0)) is the eigenvector of Q. Al- 
though the matrix Q may not be taken as a diagonal 
matrix in some convenient bases in general cases, in our 
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case the matrix Q can be diagonalized. [21] Hence, the 
diagonal elements of Q correspond to n a (x)/T directly. 
From the above discussions, all we have to do is to cal- 
culate the eigenvalues of Q. 

From the Floquet theory, the matrix Q is given by the 
eigenvalues of the Floquet Hamiltonian. [28| The Floquet 
Hamiltonian Hp is defined by 

(a', n\ Hp m) = H%~ m) + nuj5 a , (3 5 nmi (20) 

where 

H (n) = - [ He~ inut dt. (21) 
T Jo 

The \a, n) is an orthonormal basis for the matrix repre- 
sentation of Hp. The index a represents a 'state part' 
and the index n merely represents a Fourier component. 
The 'state part' means that the Hamiltonian is expressed 



as 2 x 2 matrix and hence the eigenvector of the Hamil- 
tonian has two components (states). In our case, we de- 
note the state parts as + or — (i.e., a, /3 € {+, — }). The 
| a, n) is called as a "Floquet state." Note that H is 
non-Hcrmitian, so that (a', n\ is not a conjugate state of 
| a, n), but the corresponding left orthonormal basis; we 
add the prime on a in order to clarify the fact. 



B. Perturbative calculation for the eigenvalue 

Let Hp have eigenvectors \s a<n ) and eigenvalues En.n- 
The matrix elements of Q are therefore given by [27Ll28j 

Qa0 = —£a,oS a f3- (22) 

In our case, the matrix form of the Floquet Hamiltonian 
Hp is written as follows: 



-ic\ — ic-2 + (n + 
fl(-i + l)/2 



^c_l + iC2e 



R(i 



->x 



)/2 



- ic 2 + (n + 1)uj 





R(- 
R(i 



l)/2 
-**)/2 



— ici — ic-2 + nui 
ici + ic~2e~ tx 



IC-i 






-f ic 2 e l 

- IC2 + 



(23) 



in the bases {. . . , |+, n+ 1), | — , n+ 1), +, n), — , n), . . . }. 
It is easy to see that when R = 0, we have the block di- 
agonalized matrix, and the block diagonalized elements 
have a periodic behavior. The block diagonalized Hamil- 
tonian for n = is given by 



H 



(0) 



-ici — ic-2 ic-i + ic2e 



ici + ic- 2 e 



-IC-i — IC2 



(24) 



The eigenvalues and right eigenvectors of the non- 
perturbative Hamiltonian H F °^ are denoted by e ® and 
l e a°o)' respectively. Since Hp is not Hermitian, we define 
the left eigenvectors of H^ as (s®, |. From the eigenval- 
ues of the block-diagonalized Hamiltonian H F °^ , we have 
the following total phase as the zeroth order of R: 



a4 0) (x) 



T 



(25) 



where e± x = e ±lx — 1. Equation pip indicates that the 
index + is adequate. Hence, we select jj,^ (x) as the total 
phase. It is easy to see that the total phase for R = 
case recovers the static current j (eq. ([3])) adequately. 

For i?< 1 case, the usual perturbation theory can be 
used to calculate /i+ (x) ■ There is no first order correction 



for the eigenvalues £+,o, an d hence we consider the second 
order correction: 



fc + ,0 fc 



(0) 



(26) 



where H F is the rest of the Floquet Hamiltonian Hp, 
which includes the perturbative parameter R. After the 
perturbation calculation, it is possible to calculate the 

(2) 

total current Jtotai and its fluctuation </ total as follows: 



'total 



+ $,&))' (27) 
J^ = ^{^,0(X)+^ (X)). (28) 

C. Results and discussions 

We here consider a simple case with c\ — c_i = C2 — 
c-2 = 1- In this case, the total current and its fluctuation 
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I 0.490 



0.480 



FIG. 1: (a) Total current J tota i (eq. J29}). (b) Its fluctuation 
./total ( ec l- (EOj ) . Here, we set ci = c_i = C2 = c_2 = 1, and 
i? = 0.5. 



up to the second order of R are explicitly given as 

R 2 LO 



J total 



4 1Q + uj 2 '' 



1 



R 2 



7(2) = 

total 2 16 W 



(29) 
(30) 



From these expressions, it is easy to see the effects of the 
perturbation, as follows. 

First, the perturbation induces a phenomenon like a 
stochastic resonance. Figure [IJa) shows the total cur- 
rent Jtotai calculated from eq. (f2"9"|) . Here, we set R = 0.5. 
The total current shows a peak at a certain frequency w c , 
which has been also observed experimentally. @, @] In 
the experiment in Ref. 7, Na + pumping mode of (Na,K)- 
ATPase by an oscillating electric field has been studied. 
Because different conformations of the protein have dif- 
ferent dipole moments, an oscillating electric field could 
drive structural change of the protein and hence cause 
the modulation of kinetic parameters. The actual de- 
pendency of the kinetic parameters on the external field 
could be very complicated, and we assume that two of the 
kinetic rates depend on the oscillating field in our model. 
Our model is a very simplified one, so that the quanti- 
tative comparison with the experimental results may be 



difficult. However, the qualitative behavior is the same 
as the experimental one: there is one peak, and it decays 
as ~ uj^ 1 for uj ^> 1. [9] Hence, we believe that our model 
catches the feature of the pumping phenomenon. 

Second, the perturbation varies the fluctuation of the 
total current, too. While it is possible to evaluate the 
total current by various analytical methods, [l4[ to 
our knowledge, the fluctuation has not been calculated 
explicitly for a long time. In adiabatic cases, the fluc- 
tuation has been calculated via the Berry phase inter- 
pretation. [l6[ In the adiabatic cases, it has been re- 
vealed that the perturbation decreases the fluctuation 
of the current. From eq. (f3T)]) , it is clear that the per- 
turbation decreases the fluctuation of the total current 

(2) 

even in the non-adiabatic case. The fluctuation J total is 
shown in Fig. [TJb). In contrast to the current, Fig. [T{b) 
shows no peak: the fluctuation increases monotonically 
with the frequency u, and approaches to the value of the 

( 21 

non-perturbative case (in this case, J tota i — > 1/2). 

Finally, we comment on the relationship between the 
present calculation and the Aharonov-Anandan phase. 
Because the second order correction vanishes when we 
set R = 0, one might consider that all the second or- 
der correction stems from the Aharonov-Anandan phase. 
However, it is not true. For the simplest case with 
ci = c_i = C2 = c_2 = 1, the classical current is zero, 
and the non-zero total current of eq. (|2T)|) stems from the 
Aharonov-Anandan phase. In contrast, the fluctuation 
calculated from the Aharonov-Anandan phase becomes 
zero, so that the second order term in eq. pop is induced 
from the dynamical phase. In summary, the second or- 
der correction does not stem from only the Aharonov- 
Anandan phase, but also the dynamical phase. If one 
wants to know only the 'pumping' current and its fluctu- 
ation, one more calculation for the Aharonov-Anandan 
phase is needed. [2l[ 



V. CONCLUSIONS 

In this work, we gave a unified theory for the calcula- 
tion of the total current and its fluctuation in classical 
stochastic systems under a periodic perturbation of the 
kinetic rates. The formulation is based on the full count- 
ing statistics and the non-adiabatic geometrical phase. 
It was clarified that the total current is easily obtained 
by the combination of the Floquet theory and a sim- 
ple perturbation calculation. Although it is possible to 
calculate the Aharonov-Anandan phase directly, [2l| the 
calculation is a little complicated. In a practical sense, 
the total current would be important, and as shown in 
this work, the total current is easily obtained only by the 
perturbation calculation for the eigenvalues of the Flo- 
quet Hamiltonian; the calculation is easier than the di- 
rect evaluation of the Aharonov-Anandan phase. Hence, 
the formulation given in the present paper gives a useful 
analytical method for the calculation of the current in 
the non-adiabatic cases. 
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Since our model is a simplified one, we can obtain the 
analytical expressions for the total current and its fluc- 
tuation. Although it may be difficult to compare the 
simple model and real experiments, an important theo- 
retical prediction is that there is no peak structure for 
the current fluctuation, which is different from the be- 
havior of the average current. The simple model can 
predict the qualitative behavior for the average current 
(one peak and the decay ~ w^ 1 ) adequately, so that we 
expect that the qualitative behavior for the fluctuation 
of the average current obtained in the present paper will 
be observed in real experiments. In addition, we hope 
that our analytical treatments and results become a ba- 
sis for further studies, especially for the connection with 
the non-equilibrium thermodynamics. 



where \4>(i)) is a state vector of a cyclic evolved state, 
and (4>(i)\ is related to a evolved state associated with 
the periodic adjoint Hamiltonian H'. The Aharonov- 
Anandan phase is calculated from [2(| [H| 



7(X) = ^ m)\j t m))dt. (A3) 



It has been shown that the Aharonov-Anandan phase 
gives the pump current. [2l| That is, there are the fol- 
lowing correspondences between the phases and currents: 
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APPENDIX A: FURTHER DISCUSSIONS FOR 

RELATIONSHIP BETWEEN DYNAMICAL 
PHASE AND AHARONOV-ANANDAN PHASE 

It is possible to divide the total phase fi(x) m ec l- lp"5]) 
into two parts as follows: 



(Al) 



<5(x) is called the dynamical phase, and j(x) is the 
Aharonov-Anandan phase. The dynamical phase is de- 
fined by [H 



(<t>(t)\H\<t>(t))dt, 



(A2) 



_ 1 dS( X ) 



Jc) = — 



.-/„ 



T d x 

1 dlix) 



T d x 



x=o 



(A4) 
(A5) 



The total current is given by Jtotai = */ci + */pump- Higher 
cumulants are also calculated from these phases. 

In Ref. 21, we have developed a theory for the per- 
turbative calculation for the eigenvectors of the Floquet 
state, which is related to the state vector \4>(t)). In con- 
trast, in the present paper, we clarified that the total 
current is directly related to the eigenvalues of the Flo- 
quet Hamiltonian. In practical sense, it is needed to 
evaluate the total current and its fluctuation, because 
experimental observables include these two effects. Al- 
though it is possible to evaluate the pump current from 
the Aharonov-Anandan phase, [2ll | it is more convenient 
for us to calculate the total phase directly. 
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